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Electronic analogues of nematic states, in which ro- 
tational symmetry is broken but translational invariance 
is preserved, have been proposed in a variety of correla- 
ted materials |1 1, such as quantum Hall systems |2|, cu- 
prates |3|, ruthenates |4|, heavy fermions |5| and, more 
recently, iron pnictide superconductors (H |3. In the lat- 
ter, several experiments [[6l[8l4T2t on strained samples have 
collected strong but indirect evidence that the tetragonal- 
to-orthorhombic structural transition is driven not by the 
lattice, but by electronic nematicity. However, these mea- 
surements cannot disentangled the roles of the spin lH3l - 
H6l « charge and orbital 1 17 -20J degrees of freedom in the 
nematic instability. Here, using electronic Raman spec- 
troscopy, we report the first direct measurements of inci- 
pient charge nematicity in the tetragonal phase of strain- 
free Ba(Fei_ ;E Co :E )2As2 single crystals. The diverging Ra- 
man response at low temperatures unveils an underlying 
charge nematic state that has hitherto remained unno- 
ticed. Its fluctuations persist above the superconducting 
dome, suggesting they may play a role in the pairing me- 
chanism. Comparison with the behavior of the elastic mo- 
dulus allows us to disentangle the charge contribution 
to the nematic instability, imposing strong constraints on 
theoretical models for this intricate correlated state. 

In Ba(Fei_^Co^)2As2, the structural transition at T s either 
precedes or accompanies a magnetic transition at T/v, disap- 
pearing near the doping concentration with the highest super- 
conducting transition temperature T c (see the phase diagram 
of Fig. [TJ)). The nematic/orthorhombic state is characterized 
by inequivalent Fe-Fe bond lengths along the in-plane a and b 
directions (x and y coordinates respectively of the one-Fe unit 
cell used throughout, see Fig.[T^), and by anisotropic electro- 
nic responses lETTl . If this state is a consequence of the nema- 
ticity of electronic degrees of freedom, as is widely held, as- 
sociated fluctuations should be present in the tetragonal phase 
and should increase as the temperature is lowered towards T s . 
Probing these nematic fluctuations in the magnetic channel 
is a difficult task as it involves measuring higher order spin 
response functions. On the other hand, the charge nematic 
fluctuations, if present, can be extracted via electronic Raman 
scattering measurements. 

The electronic Raman response probes the weighted charge 
correlation function (p^(uj)p^(-uj)}, where p M = Xlk7k n k 



depends on the charge-density operator of the momentum 
state k, and on the form factor 7^ whose symmetry /i is de- 
termined by the polarizations e/ and of the incident and 
scattered photons Q. To probe the in-plane charge nematic 
fluctuations, two polarization configurations can be conside- 
red (see inset of Fig.[ljl). For photons polarized along the dia- 
gonals of the Fe-Fe bonds, the form factor has x 2 — y 2 (Bi g ) 
symmetry, and is sensitive to nematic order along the Fe-Fe 
bonds. This is the type of C4 (tetragonal) symmetry -breaking 
realized in the iron pnictides. For photon polarizations paral- 
lel to the Fe-Fe bonds, the form factor has xy (i?2#) symme- 
try, being sensitive to nematic order along the diagonal of the 
Fe-Fe bonds. Although this kind of symmetry breaking is not 
observed in the iron pnictides, it has been suggested to oc- 
cur in URu2Si2, which has the same crystalline structure as 
BaFe2As2 [ 5 ]. The behaviors of these form factors in momen- 
tum space are depicted in Fig.[lJ. 

Due to the fluctuation-dissipation theorem, the dynamic 
charge nematic fluctuations should be manifested in the ima- 
ginary part of the Raman response function (x^)", which is 
directly measured in Raman scattering experiments (see Sup- 
plementary Information), in the appropriate symmetry /i, na- 
mely, the x 2 — y 2 symmetry l23l l24l . This is illustrated in 
Fig. [TJi for a twinned, stress-free, single crystal of the parent 
compound BaFe2As2, where (x M )" is plotted as function of 
frequency for different temperatures and for the two symme- 
tries described above. While the response in the xy symmetry 
is essentially temperature independent above T s = 138 K, the 
x 2 — y 2 response displays a considerable build-up of inten- 
sity below 500 cm -1 upon approaching T s , with a subsequent 
collapse in the nematic/orthorhombic phase. The temperature 
dependence and the distinctive x 2 — y 2 symmetry of this low 
frequency quasi-elastic peak (QEP) clearly links it to dynamic 
charge nematic fluctuations corresponding to an orientational 
order along the Fe-Fe bonds. 

While the spectral lineshape of the QEP is linked to the 
relaxational dynamics of the nematic fluctuations, we choose 
here to concentrate on a more transparent quantity : the sta- 
tic charge nematic susceptibility. Indeed the strong increase 

of the QEP intensity can be linked to a diverging static charge 
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nematic susceptibility, Xo ~ V » Nm me Kramers-Kronig rela- 
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Figure 1: Raman scattering measurements of charge nematic fluctuations in BaFe2As2 a, Tetragonal FeAs layer, with the x and y 
axes defined along the Fe-Fe bonds b, Sketch of the phase diagram of Ba(Fei_ x Co x )2As2. T s , Tat and T c are the structural, magnetic and 
superconducting (SC) transition temperatures respectively, c, Momentum-space structure of the form factor 7^ for x 2 — y 2 and xy symmetries 

(see Supplementary Information), d, Temperature dependent Raman response (\ x ~ v )" and (x xy )" m a twinned BaFe2As2 single crystal 
with T s — 138 K. The incoming and outgoing photon polarizations (ei ,es) used for each symmetry configuration are depicted in the insets. 
The sharp peaks are due to phonon excitations. The electronic Raman continuum in x 2 — y 2 symmetry displays a low frequency quasi-elastic 
peak (QEP) that is superimposed on a weaker and broad continuum that extends to energies above 1000 cm -1 and is essentially temperature 
independent in the tetragonal phase. In the orthorhombic phase, this broad continuum shows a suppression below 500 cm -1 in both symmetries 
because of the Fermi surface reconstruction induced by the simultaneous magnetic order [25] . 



tion (see Supplementary Information) : 

O poo 

Xo~ y =~ duj( X "f- y2 (u)/u (1) 

The relevant quantity governing the static nematic suscep- 
tibility is thus the Raman conductivity x" highlighting 
the importance of the low frequency part of x" m determi- 
ning Xo ~ V • The temperature dependence of x" '/w, where 
the QEP is now centered at zero frequency, is shown in Fig. 
[5^ for six different Co concentrations of Ba(Fei_ :c Co :E )2As2, 
spanning the phase diagram from the parent x = composi- 
tion (T s = 138 K and T c = 0) up to the strongly overdoped 
x = 0.20 composition (T a = T c = 0). For x < 0.045, the 
QEP displays a systematic enhancement as temperature is lo- 
wered towards T s before collapsing in the symmetry broken 
phase. The diverging behavior of the QEP extends down to T c 
for x = 0.065 where the superconducting transition tempera- 
ture is optimal and no structural transition is detected. For this 
particular composition the QEP was found to disappear qui- 
ckly upon entering the superconducting state indicating a sup- 



pression of nematic fluctuations in the superconducting state 
(not shown). Above optimal composition, the enhancement 
of the QEP is strongly reduced but remains sizable even for 
x = 0.10, before disappearing for x = 0.20 where x" i s 
essentially a flat, temperature independent background. 
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The static nematic charge susceptibility Xo v was extrac- 
ted using equation [9] via a partial integration of the Raman 
conductivity up to 500 cm -1 , since above this frequency the 
spectra are temperature independent in the tetragonal phase. In 
order to perform the integration the Raman conductivity spec- 
tra were interpolated linearly from the lowest frequency expe- 
rimentally accessible, 9 cm -1 , down to zero (Supplementary 
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Fig. 1). The doping and temperature dependence of Xo V 
are summarized in the phase diagram of Fig. [5]). The maxi- 
mum of the static charge nematic susceptibility closely tracks 
the structural transition temperature in the underdoped region, 
vanishing near optimal doping. This temperature and doping 
dependence is qualitatively consistent with previous anistro- 
pic transport data of strained detwinned crystals |6, 8] [TOl. 
However, resistivity anisotropy is only an indirect probe of 
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Figure 2: Doping dependence of charge nematic fluctuations and static charge nematic susceptibility in Ba(Fei_ x Co x )2As2 a, Tempe- 
rature dependent Raman conductivity (\ x ~ v )" /oj for x = (parent), x = 0.02 (slightly underdoped), x = 0.045 (underdoped), x = 0.065 
(optimally doped), x = 0.10 (overdoped), and x = 0.20 (strongly overdoped). The structural transition temperature, as determined in-situ both 
by looking at phonon anomalies [25 1 and by the onset of Rayleigh scattering by structural twin domains, is indicated for the three underdoped 
compositions. The x = 0.065 composition corresponds to optimal superconducting transition temperature (T c = 24.5 K) where no structural 
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transition was detected, b, Evolution of the static charge nematic susceptibility, Xo y , as a function of temperature and doping. The struc- 
tural transition temperature T s and the superconducting transition temperature T c are indicated in red squares and blue triangles respectively. 
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Xo y displays a clear maximum at the structural transition temperature. 



the nematic order parameter since it cannot disentangle the 
various possible sources of electron nematicity, and, further- 
more, is sensitive to scattering mechanisms l2T1l . By contrast 
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Xo ~ v is a direct measure of charge nematic fluctuations and 
is non-zero even in strain-free samples where the C4 symme- 
try is preserved. 

To perform a more quantitative analysis, in Fig. [3^ we 
plot the inverse susceptibility as a function of temperature 
in the tetragonal phase (T>T S ) and for the six Co composi- 
tions. The diverging tendency, manifested as a softening of 
the inverse susceptibility, is seen for all compositions up to 
x = 0.10, being absent only for the strongly overdoped, non- 
superconducting, x = 0.20 composition. For all other compo- 
sitions the inverse susceptibility above T s can be well descri- 
bed over a large temperature range, spanning at least 150 K, 
by a simple Curie-Weiss law of the form : 

( x f-yy\T) = C-'^f (2) 
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where C is the high temperature limit of Xo V > ^0 is the 
mean-field charge nematic transition temperature, and 6 < T 
is a constant. The resulting fits for the inverse susceptibility 
are shown in Fig. [3^. They unveil an incipient charge ne- 
matic instability at To over a wide doping range, which in- 
cludes the superconducting dome, in the phase diagram of 
Ba(Fei_ a ,Co :c )2As2. 



The extracted To follows the trend of the thermodynamic 
structural transition temperature T s , decreasing with doping 
and vanishing near x ~ 0.06. However, T is significantly 
smaller than T s , by about 50 K, across the entire phase dia- 
gram (see Fig. [3])). This fact, in conjunction with the obser- 
vation of the build-up of the charge nematic fluctuations over 
a large temperature range, allow us to conclude that the inci- 
pient charge nematicity is not a mere consequence of a cou- 
pling with the lattice orthorhombicity. More importantly, since 
the Curie- Weiss expression ^ with To significantly smaller 
than T s describes very well the data up to a few Kelvin above 
T 39 it implies that the incipient charge nematicty is, in fact, 
weakly coupled to the lattice. Here, we emphasize that tetra- 
gonal symmetry-breaking has to occur at the same tempera- 
ture in both elastic and charge degrees of freedom. Therefore 
our analysis also suggests the presence of another nematic de- 
gree of freedom which drives the structural transition at T s 
and preempts the charge nematic transition at a lower tempe- 
rature To . 

These conclusions can be reached in a more transparent 
way, without relying on any particular fitting expressions, sim- 
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ply by comparing the Xo V data directly with the shear mo- 
dulus C s = C\\ — C12, which is the orthorhombic lattice 
stiffness ElEHl. In fact, in Ref. CD a Curie-Weiss expres- 
sion analogous to Eq. (pi was used to fit the softening of the 

(c ) 

shear modulus, finding Tq sJ w T s . Because symmetry re- 
quires the charge nematic order parameter to couple linearly 
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Figure 3: Curie- Weiss analysis of the inverse static charge ne- 
matic susceptibility in Ba(Fei_ cc Co cc ) 2 As2 a, Temperature depen- 
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dence of the inverse nematic charge susceptibility, (xo V )~ » m 
the tetragonal phase (T>T S ) as a function of Co composition. The 
lines are Curie- Weiss fits for each composition (see text), b, (x, T) 
phase diagram showing the orthorhombic (Ortho) and superconduc- 
ting (SC) phases. The mean-field transition temperature extracted 
from the Curie-Weiss fit, To, is shown in green square (the green line 
is a linear fit of its doping dependence). The corresponding structural 
transition temperature T s (red squares), magnetic transition tempera- 
ture T/v (white circles), and superconducting transition temperature 
T c (blue triangles) are also indicated |26|. 



to the orthorhombic distortion, there is a simple relationship 
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between C s and Xo ~ V if charge nematicity is the only soft 
mode present H6]|23 : 



eg 



1 



eg 



x — y 

Xo 



(3) 



Here, A is the coupling to the lattice and C^ is the high- 
temperature shear modulus. In Fig. [4] we test the validity 
of the above relation for the parent compound BaFe2As2 by 
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Figure 4: Charge nematicity and lattice softening in BaFe2As2. 

Temperature dependence of the experimental shear modulus (data 
of Ref. (27), in purple), together with the expected temperature de- 
pendence of the shear modulus due to the coupling between the the 
charge nematic and orthorhombic order parameters using equation [3] 
(in black). The bare shear modulus was assumed to be temperature 

A 2 

independent and the only adjustable parameter, a = (see text) 

was chosen so as to fit the shear modulus data at high temperatures 
(a — 0.8). The temperature scale was normalized using the measu- 
red structural transition temperatures. Notice that the lattice softening 
induced by the charge-nematic fluctuations is insufficient to describe 
the measured softening, implying the existence of additional electro- 
nic soft modes. 



comparing the C s values inferred from our Xo ~ y data vm 
equation [3] (black) with the experimental C s data of Ref. l27l 
(purple). It is clear that the relationship ([3} cannot be satisfied, 
i.e. the observed lattice softening cannot be accounted solely 
by a linear coupling with the charge nematic order parameter. 
Consequently, there must be an additional electronic nematic 
degree of freedom that not only drives the structural transition, 
but also induces the charge nematic softening. 

One possibility is that spin is this additional degree of free- 
dom. Indeed, theoretical models considering either localized 
[Q31IT4) or itinerant spins E5l[T6 | have proposed that magne- 
tic fluctuations lead to a spontaneous electronic C4 symmetry 
breaking already in the paramagnetic phase. Alternatively, it 
has also been proposed that the structural transition is driven 
by orbital ordering between xz and yz 3d Fe orbitals fTTvl— 
l20l . Although the charge fluctuations measured here do not 
necessarily come only from fluctuations of the relative charge 
nxz — nyz between these two orbitals, it follows from the or- 
bital content of the Fermi surface of the iron pnictides that 
these orbital fluctuations should give a major contribution to 

2_ 2 

Xo y if orbital order is the driving instability (|2|, see also 
Supplementary Informations). 

Thus, our observation, that charge nematic fluctuations 
alone cannot drive the structural transition, imposes strong 
constraints on models of nematicity in iron-pnictides. The 
persistence of these fluctuations above the entire supercon- 
ducting dome raises the question of whether they play a role 
in the pairing mechanism l2T1l . Besides shedding light on 
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the nature of the nematic state of the pnictides, our approach 
provides a novel route to investigate electronic nematicity in 
other strongly correlated systems where this type of state has 
been proposed, such as the magnetic field induced nematic 
phase of Sr 2 Ru307 |4 |, the hidden-order phase of URu2Si2 
l5l . and the pseudogap in underdoped cuprates l30l . 
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Supplementary information 

A. Methods 

Raman experiments have been carried out using a diode- 
pumped solid state laser emitting at 532nm and a triple grating 
spectrometer equipped with a nitrogen cooled CCD camera. 
Special care was taken in order to determine laser heating. It 
was first estimated by comparing the power and temperature 
dependences of the phonon frequencies. This estimate was 
then cross-checked by monitoring the onset of Rayleigh scat- 
tering by orthorhombic structural domains across the structu- 
ral transition temperature as a function of laser power. Both 
methods yielded an estimated heating of 1 K ± 0.2 per mW 
of incident power. In order to extract the imaginary part of 
the Raman response function, the raw spectra were corrected 
for the Bose factor and the instrumental spectral response (see 
also Supplementary information B). 

Single crystals of Ba(Fei_ :E Co a ;)2As2 were grown using 
the self-flux method. The Co content was checked using an 
electron probe (Camebax 50) yielding a Co content within 
0.5%. The magnetic transition temperature was determined by 
transport measurements performed on crystals from the same 
batch [26]. The superconducting transition temperature of the 
crystals were measured by SQUID magnetometry. The struc- 
tural transition temperature was determined by monitoring the 
onset of Rayleigh scattering by structural twin domains at very 
small incident laser power. 

B. Extraction of the static charge nematic susceptibility from 
Raman scattering measurements 

The experimentally measured Raman intensity in the sym- 
metry fx, is proportional to the weighted charge corre- 
lation function 

I^(uj) oc S»(u) = (p^(uj)p^(-uj)} (4) 

The correlation function is in turn directly linked to the 
imaginary part of the Raman response (x^)" via the fluctua- 
tion dissipation theorem 



S"M = — (l + n( W ,T))(x")"H (5) 

7T 

where n(u,T) is the Bose-Einstein distribution function. The 
electronic Raman response function x M is given by 

X»H = -j^ dte»*(\pl>(t),ff(0)]) (6) 
where the operator p M has the form 

^ = £ 7 £n k (7) 

k 



in terms of the charge density operator in the momen- 
tum state k, and a form factor, also called the Raman vertex, 
7k( e i? e s)> whose symmetry index fi depends on the polari- 
zations of the incident and scattered lights and e s respec- 
tively. 1 1 ] Since the photon wave vector is several orders of 
magnitude smaller than the typical Brillouin zone size, there 
is negligible momentum transfer in the electron-photon scat- 
terings, and therefore Raman spectroscopy probes the system 
uniformly. 

The x 2 — y 2 or B\ g symmetry can be selected by choosing 
crossed incoming and outgoing photon polarizations at 45 de- 
grees with respect to the Fe-Fe bonds. Here the notation B\ g 
refers to the one Fe unit cell whose axes are along the Fe-Fe 
bonds. Similarly the xy or B 2g symmetry can be selected by 
choosing incoming and outgoing photon polarizations along 
the Fe-Fe bonds. Note that in terms of the full lattice unit cell 
(or 2 Fe unit cell), which has its axes at 45 degrees to the Fe- 
Fe bonds and is more commonly used in the Raman literature, 
the Big (B2 9 ) symmetry discussed here corresponds to the 
B 2g (Big) symmetry. 

The momentum space structure of 7^ is constrained by 
symmetry. For example in the case of x 2 — y 2 symmetry, 7^ 
must change sign under mirror symmetry with respect to the 
direction at 45 degrees of the x and y axis. Using the ef- 
fective mass approximation for a tight binding model with 
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nearest neighbour hopping integrals only we have 7^ y = 
cosk x — cosky and 7^ = sink x — sink y dill. 

The k- space structure of these form factors are shown in 
Fig. 1 of the main manuscript. While we only have access to 
the the imaginary part of the symmetry dependent response 
as a function of frequency, we can extract the corresponding 
static susceptibility, Xo , using Kramer-Kronig relation linking 
the real and the imaginary parts of the Raman response func- 
tion : 




Taking the u = and using the fact that x" is an °dd func- 
tion of oj, we obtain the following expression for the symmetry 
dependent static susceptibility : 

Xff = - / duj\ X »)"(uj')/uj' (9) 

The expression above shows that the relevant quantity 
controlling the static susceptibility is not the Raman response 
x" but the Raman "conductivity" x" which is dominated 
by its low frequency behavior. 

Experimentally only the response in x 2 — y 2 symmetry, 
which measures charge nematic fluctuations, shows a signi- 
ficant build-up at low frequency or quasi-elastic peak (QEP) 
upon cooling. The QEP appears as a peak centered at zero- 
frequency in the raw intensity data, /, while it is pushed to 
finite frequency when the quantity x" is plotted as in the main 
manuscript (see equation [5} Since we have only access to a 
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FIGURE 5: Linear extrapolation of (x x ~ v )" /w fox x = (pa- 
rent compound BaFe2As2) and for 10 different temperatures in the 
tetragonal phase. The linear extrapolation was performed using the 
frequency window 9-25 cm -1 . 



finite frequency range, the integral in Eq. [9] can only be per- 
formed up to a finite frequency cut-off. Experimentally the 
QEP above T s is superimposed on a weaker and broad electro- 
nic continuum that extends up to energies above 1000 cm -1 . 
This continuum shows a reconstruction in the SDW state (4) 
but is essentially temperature independent in the tetragonal 
phase. In particular \" I s t0 within our experimental accu- 
racy temperature independent above 500 cm -1 in the tetrago- 
nal phase. One can therefore reliably extract the temperature 
dependent charge nematic susceptibility by restricting the in- 
tegral to energies lower than 500 cm -1 where the response 
in temperature dependent and dominated by the QEP. On the 
low energy side, the Raman measurements were performed 
down to 9 cm -1 . In order to perform the integration down 
to zero frequency, the Raman conductivity, (\ x ~ y )"{uj')/uj' 
was extrapolated linearly down to uj = as shown in Fig. [5] 
The resulting static nematic susceptibility is shown as a func- 
tion of temperature and doping in Fig. 2 of the manuscript. 



C. Raman Vertex in the orbital basis 

In a multi-orbital system, the basic operator for Raman res- 
ponse in the symmetry channel fi is given by : 



(10) 



kcr m,n 



where kcr creates an electron with momentum k and spin 
a in orbital m. For the x 2 — y 2 channel, ignoring screening 
effects and vertex corrections, the free-electron form factor or 
Raman vertex depends only on the band dispersion e mn (k) : 



l X mn V (k) 



d 2 r 



(k) 



3 2 f 

U Con 



(k) 



dkl 



dkl 



(11) 
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FIGURE 6: x 2 — y 2 form factor (cos k x — cos k y ) in grey scale. A 
sketch of the typical Fermi surface sheets of the iron pnictides has 
been superimposed |3 1. Their orbital content is also indicated. 



In the iron pnictides, all 3d orbitals may contribute to the 
Fermi surface. Let us focus on the role played by the orbitals 
xz and yz. In the tetragonal phase, their intra-orbital disper- 
sions are identical upon a 90° rotation of the coordinate sys- 
tem, i.e. 



&xz,xz (foxi hy) — &yz,yz ( ^yi ^x) (12) 

Therefore, it follows for the intra-orbital Raman vertices : 



Ixz.xz (k x ,k y ) — Jyz.yz ( k y ,k x ) 



(13) 



At low energies, the main contribution to the Raman scatte- 
ring comes from electronic states at the Fermi level. Symme- 
try requires that if a point (k x , k y ) at the Fermi surface has xz 
orbital character, then the point at (—k yi k x ) also belongs to 
the Fermi surface and has yz orbital character. Combined with 



the symmetry relation ( 13 ), this implies that the charge den- 
sity difference Sn = n xz — n yz , with the proper overall form 
factor, contributes to the effective x 2 — y 2 Raman charge. 

Similar conclusions can be drawn by analyzing the orbi- 
tal content of the Fermi surface of the iron pnictides and the 
x 2 — y 2 form factor cos k x — cos k y . In Fig. [6j we superim- 
pose a sketch of the typical Fermi surface of the iron pnictides 
to the form factor plotted in Fig. lb of the main text. It is 
clear that the form factor changes sign between the electron 
pockets located at X = (ir, 0) and Y = (0, 7r). First principle 
calculations Q reveal that while the X electron pocket has 
mostly d yz character (and no d xz contribution), the Y pocket 
has mostly d xz character (and no d yz contribution). Therefore, 
since the x 2 — y 2 form factor changes sign between the X and 
Y pockets, and they have symmetry-related d xz and d yz spec- 



tral weights, the relative charge Sn between the d xz and d, 
orbitals appears in the x 2 — y 2 response. 



yz 
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Notice that the remaining 3d orbitals, as well as inter-orbital 
terms, in principle, also contribute to the x 2 — y 2 Raman res- 
ponse. A detailed discussion of all the non-zero intra-orbital 
and inter-orbital Raman vertices 7^ n (k) was presented in 
Ref. |2|. For instance, the xy intra-orbital vertex also satisfies 
the relationship : 

lxy,xy (kx-> k y ) = — ^xy^xy (~^y^x) (14) 



implying that a breaking of tetragonal symmetry could in prin- 
ciple be driven solely by the xy orbital. However, if the ne- 
matic instability is driven by a spontaneous d xz /d yz orbital 
polarization, as suggested by a few theoretical models, then 
the fluctuations associated with the charge density difference 
8n between xz and yz orbitals (Sn 5n) should dominate the 
x 2 — y 2 Raman susceptibility. 
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